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ABSTRACT 

We report on recent results on the decay of the (j) into tt+tt^ and radiative 
decays into 7r°7r''7 and ir^rj^, which require the consideration of the final state 
interaction of a pair of mesons in a region inaccessible to Chiral Perturbation 
Theory. By using nonperturbative chiral unitary methods for the meson meson 
interaction we can obtain the corresponding decay widths and the results are 
compared with recent experimental data. 



1 The 4> — > TT+TT decay 

The (j) decay into vr+vr^ is an example of isospin violation. The (j) has isospin 
/ = 0, spin J = 1, and hence it does not couple to the tt+tt^ system in the 
isospin limit, which implies the rule I + J — even. The experimental situation 
on this decay is rather confusing. There are two older results whose central 



values are very different but their quoted errors are so big that both were still 
compatible: The first one from 11) gives BR = (1.94+1.03-0.81) x 10"*. The 
second one from i) provides BR = (0.63 + 0.37 - 0.28) x 10""*. Very recently 
two new, more precise, but conflicting results have been reported from the two 
experiments at the VEPP-2M in Novosibirsk: the CMD-2 Collaboration reports 
a value BR = (2.20 ± 0.25 ± 0.20) x lO"'^ & whereas the SND Collaboration 
|) obtains BR = (0.71 ± 0.11 ± 0.09) x lO^-^. 

Isospin violation has become a fashionable topic in Chiral Perturbation 
Theory (xPT) but the ^ TTTT decay is however unreachable with plain 

X-PT, since it involves the propagation of the pair of pions around 1 GeV, far 
away from the xPT applicability range. 

Nevertheless, new nonperturbative schemes imposing unitarity and still 
using the chiral Lagrangians have emerged enlarging the convergence of the 
chiral expansion. In 0) the inverse amplitude method (lAM) is used in one 
channel and good results are obtained for the cr, p and K* regions, amongst 
others, in tttt and ttK scattering. In H the method is generalized to include 
coupled channels and one is able to describe very well the meson-meson scat- 
tering and all the associated resonances up to about 1.2 GeV. A more general 
approach is used in 0) by means of the N/D method, in order to include the 
exchange of some preexisting resonances explicitly, which are then responsible 
for the values of the parameters of the fourth order chiral Lagrangian. 

Here we shall follow the work H) since it provides the most complete 
study of the different meson-meson scattering channels, including the mesonic 
resonances and their properties up to 1.2 GeV. In particular, this method yields 
a resonance in the / = 0, J = 1 channel, the ws resonance, related to the 4>, 
and this allows us to obtain the strong contribution to the — > tttt decay. We 
also consider electromagnetic contributions at tree level which turn out to be 
dominant and were already considered in 

In order to calculate the contribution of an intermediate photon to the 
(j) — > TTTT decay, let us consider the effective Lagrangian for vector mesons pre- 
sented in lij' , which is written in terms of the SU(3) pseudoscalar meson 
matrix $ and the antisymmetric vector tensor field V^i^ defined in 



c,[v{i- -)] = ^{v,,.nn + '-^{v.^u^un, (1) 



where "( )" indicates the SU(3) trace. In order to introduce the physical states 



(j) and uj, we assume ideal mixing between the loi and vector resonances and 
hence taking into account that the uji does not couple to pairs of mesons at 
the order of eq. (1), the coupling of the (f) is easily deduced from that of the 
cjs by simply multiplying the results of the by the factor — . With these 
ingredients and the standard 77r7r coupling we can write the contribution of a 
Feynman diagram with the 4> going to a photon which then couples to a pair 
of pions, and which is given by 

iC^^^^-=ie^^^^e>^{4>){p+ ^p^)^F{Ml), (2) 

where p+ and p- are, respectively, the momenta of positive and negative pions 
and F{q^) is the pion electromagnetic form factor, which at the (j) mass is given 
by i^(M|) = —1.56 + i 0.66. This can be compared with the coupling of the 
to if + , or K'^ K'^, which can be obtained from the Gy term in Eq. (|l|) and 
reads 

i^^K+K- = -ig^K+K- e-^{4>) {P+ -P-)ti , 9^K+K- = ^Ij^p ' 

which provides the right (f) decay width with a value of Gy = 54.3 MeV . 
By analogy to Eq. (||), Eq. (||) gives a 4> couphng to tt+tt" 



/)7r+7r 



which provides the tt+tt" decay width with the tree level photon mecha- 
nism. With a value of Fy — 154 MeV from the p e+e" decay 111) and using 
the coupling of Eq. (^) one obtains a branching ratio to the total width of 
1.7 X 10^"*. 

In order to evaluate the strong contribution to the process we consider the 
KK — > TT+TT" amplitude corrected from isospin violation effects due to quark 
mass differences. The method used is based on the chiral unitary approach 
to the meson-meson interaction followed in H & . The technique starts from 
the 0{p^) and O(p^) X-PT Lagrangian and uses the lAM in coupled channels, 



ouph 

n. 



generalizing the one channel version of the lAM developed in 

Within the coupled channel formalism, the partial wave amplitude is given 
in the lAM by the matrix equation 



r-T2[r2-T4]-iT2, 



(5) 



where T2 and T4 are Oip^) and 0{p^) xPT partial waves, respectively. In 
principle T4 would require a full one-loop calculation, but it was shown in H) 
that it can be very well approximated by 



where is the tree level polynomial contribution coming from the £4 chiral 
Lagrangian and G is a diagonal matrix for the loopfunction of the intermediate 
two meson propagators which are regularized in Ef by means of a momentum 
cut-off. 

In the present case, in which isospin is broken explicitly and J — 1, we are 
dealing with three two-meson states: K+R-, K°K° and tt+tt", that we will 
call 1, 2 and 3, respectively. The amplitude is a 3 x 3 matrix whose elements 
are denoted as Ty. The T2 and T^" amplitudes used in the present work and 
calculated in the isospin breaking case, are collected in the appendix of 0). 
The fit of the phase shifts and inelasticities is carried out here in the isospin 
limit, as done in H) and there are several sets of Li coefficients which give rise 
to equally acceptable fits. 

We write in table |l| the values of the coefficients of the different sets of chi- 
ral parameters. The corresponding results for the phase shifts and inelasticities 
can be seen in 0) where it is shown that the small differences in the results 
appear basically only in the ao(980) and k(900) resonance regions, where data 
have also larger errors or are very scarce. 

In order to evaluate the contribution to the cj) tt+tt^ coupling from the 
strongly interacting sector we evaluate the K^K^ K^K^ amplitude (Tn) 
and the K^K~ tt+tt" amplitude (Tia) near the pole of the cos resonance 
which in our case appears around M^^ = 920 MeV. Close to the pole the 
amplitudes obtained numerically are then driven by the exchange of an wg- 

By assuming a coupling of the type of Eq. (^) for the ujg to K^K^ and 
TT+TT^, these two amplitudes, close to the pole, are given by 



Rer4 ~ tI +T2ReGT2 



(6) 



Til 



9, 



'ujsK+K- p2 _ yT2 



^VK -PK' 



1 



^PK -Ptt- 



(7) 



T13 — 9u}sK+K- 5w87r+7r 




where pi is the three-momentum of the i particle in the CM frame. 





Li 


L2 


La 


L4 


Ls 


2L6 + ^8 


L7 


Qmax 




set 1 


0.91 


1.61 


-3.65 


-0.25 


1.07 


0.58 


-0.4 


666 MeV 


1.3x10--* 


set 2 


0.91 


1.61 


-3.65 


-0.25 


1.07 


0.58 


0.05 


751 MoV 


1.0x10-* 


set 3 


0.88 


1.54 


-3.66 


-0.27 


1.09 


0.68 


0.10 


673 MeV 


1.3x10"* 




Li 


L2 


is 


Li 


is 


2Le + Ls 


Lr 






ChPT 


0.4 
±0.3 


1.4 
±0.3 


-3.5 
±1.1 


-0.3 
±0.5 


1.4 
±0.5 


0.5 
±0.3 


-0.4 
±0.2 




ref. 0) 





Table 1: Different sets of chiral parameters (xlO^^ ) that yield reasonable 
fits to the meson-meson scattering phase shifts and the corresponding cj) —t tttt 
branching ratio prediction. We have used a hat to differentiate them from those 
obtained for standard ChPT. However, as it is explained in II) , we still expect 
them to be relatively similar once the appropriate scales are chosen (roughly 
ji ~ \.2qmax, see I3i for details). 



By looking at the residues of the amphtudes Tn, T13 in the pole we 
can get the products g^K+K- 9,pK+K- and g^x+K- 9<i>^+7r-- Thus, defining 

Q.u= lim (p2_M2)-i^ (8) 

we obtain the ratio of the g,f,K+K- to g^x+K- by means of the ratio of Qi^ 
to Qii, and hence taking g^K+K- from Eq. we get the value for g'"^^+^_. 
Then, by adding the above contribution with that of Eq. (||) we can obtain the 
^ TT+Tf- decay width. We have taken Fy Gy > 0, as demanded by vector 
meson dominance EH^. 

Each set of chiral parameters has then been used in the isospin-breaking 
amplitudes given in the appendix of , obtaining a value of BR{(j) tttt) 
given in table The dispersion of the results provides an estimate of the 
systematic theoretical uncertainties. 

From table |l|, we obtain, after taking into account the strong contributions 

Br{(f) — > 7r7r)trce-|-strong — (1.2±0.2) X 10^'' (9) 

On the other hand, explicit calculations of the absorptive part of the 777 in- 
termediate channel 0) give a contribution of about 1/4 of the kaon loops. 
In order to estimate the uncertainties from neglecting the photonic loops we 
take a conservative estimate and consider them of the same magnitude as the 
strong interaction correction, and, hence, add an extra ±0.5 x 10"^ uncertainty. 




(d) (e) 



Figure 1: Diagrams for the decay (f> — s- tt^tt^j. 



Adding in quadrature the errors from the different sources, our final result is 
the band of values: 

^ tttt) ~ 0.7to f.7 X 10"'', (10) 

which is compatible with the present PDG average within errors and lies just 
between the results of the two recent experiments, which are much more precise, 
but mutually incompatible. 

2 The radiative decay into 7r°7r°7 and Tr^rj-f 

The (p meson cannot decay into two pions or tt'^t] in the isospin limit. The 
decay into two neutral pions is more strictly forbidden by symmetry and the 
identity of the two pions. As a consequence the decay of the (p into 7r°7r°7 and 
TT^rjj is forbidden at tree level. However, the ip decays into two kaons and the 
processes described can proceed via the loop diagrams depicted in Fig. 1 where 
the intermediate states in the loops stand for KK. 

The evaluation of the diagrams of Fig. 1 is done in . The terms with 
Gv of Eq. (2) contribute to all the diagrams in the figure. However, the Fy 
term of Eq. (2) only contributes to the diagrams containing the contact vertex 
(j) jKK, like diagrams (a), (e). The idea follows closely the work of 
but for the treatment of the final state interaction of the mesons one uses here 
the norperturbative chiral techniques. In this case for L=0, which is the only 
partial wave needed, one can use the results of Elt , where it is proved that the 



use of the Bethe Salpeter equation in connection with the lowest order chiral 
Lagrangian and a suitable cut off in the loops gave a good description of the 
meson meson scalar sector. Furthermore, in it was proved that the meson 
meson amplitude in those diagrams factorized on shell. The loops of type (a), 
(b) and (c) can be summed up using arguments of gauge invariance following 
the techniques of and lead to a finite amplitude. On the other hand, 

the terms involving Fy and a remnant momentum dependent term from the 
Gy Lagrangian in Eq. (2) only appear in the contact vertex jKK^ and 
the diagrams of type (b), (c) are now not present. Hence, in this case the only 
loop function involved is the one of two mesons which is regularized as in 
for the problem of the meson meson scattering. The average over polarization 
of the (j) for the modulus square of t matrix is then easily written and for the 
case of 7r''7r°7 decay one finds 



For the 



PV3 
TT^rij case we have 
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where Gr+k- and Gk+k~ are the loop functions mentioned above. 

We have evaluated the invariant mass distribution for these decay chan- 
nels and in Fig. 2 we plot the distribution dB/dMj for — > 7r°7r*'7 which 
allows us to see the — > /07 contribution since the /o is the important scalar 
resonance appearing in the K^K^ tt'^tt" amplitude EH^. The results are 
obtained using G'y=55 MeV and Fv/ = 165 MeV, which are suited to describe 
the KK and e+e~ decay of the 0. The solid curve shows our prediction, with 
FyGv > 0, the sign predicted by vector meson dominance, as we quoted above. 
The dashed curve is obtained considering FyGy < 0. In addition we show also 
the results of the intermediate dot-dashed curve which correspond to taking 
for Gv and Fy the parameters of the p decay, Gv=69 MeV and fy = 154 MeV. 
We compare our results with the recent ones of the Novosibirsk experiment 
HI) . We can see that the shape of the spectrum is relatively well reproduced 
considering statistical and systematic errors (the latter ones not shown in the 
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Figure 2: Distribution dB / dMi for the decay — > with Ml the invariant 

mass of the system. Solid line: our prediction, with FyGy > 0. Dashed 

line: result taking FyGy < 0. The data points are from El^ and only statistical 
errors are shown. The systematic errors are similar to the statistical ones . 
The intermediate, dot- dashed curve corresponds to the results obtained using the 
Gv and Fy parameters of the p decay. 



figure). The results considering FyGy < are in complete disagreement with 
the data. 

The finite total branching ratio which we find for the tt'^tt'^j decay is 
0.8 X 10^^, which is slightly smaller than the result given in 111, (1.14±0.10± 
0.12) X lO"'', where the first error is statistical and the second one systematic. 
The result given in H) is (1.08 ± 0.17 ± 0.09) x 10^**, compatible with our 
prediction. Should we use the values for Fy and Gy of the p decay we would 
obtain 1.7 x 10~*. The branching ratio obtained for the case cj) n'^rjj is 
0.87 X 10""*. The resuhs obtained at Novosibirsk are H) (0.83 ± 0.23) x lO"'^ 
and 11) (0.90 ± 0.24 ± 0.10) x lO"'^. Should we use the values for Fy and 
Gy of the p decay we would obtain 1.6 x 10^^. The spectrum, not shown, is 
dominated by the oq contribution. 

The results reported here are two examples of the successful application 



of the chiral unitary techniques. A recent review of multiple applications of 
these methods can be seen in H) . 
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